A Fermi Fluid Description of the Half-Filled Landau Level 
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We present a many-body approach to calculate the ground state properties of a system of electrons 
in a half-filled Landau level. Our starting point is a simplified version of the recently proposed 
trial wave function where one includes the antisymmetrization operator to the bosonic Laughlin 
state. Using the classical plasma analogy, we calculate the pair-correlation function, the static 
structure function and the ground state energy in the thermodynamic limit. These results are in 
good agreement with the expected behavior at ^ = i . 



The fractional quantum Hall effect (FQHE) which 
is understood to be due to condensation of electrons to 
unique incompressible states as a result of electron cor- 
relations j^,^, fails to explain the odd behavior of even- 
denominator filling fractions which lie right in the mid- 
dle of all the observed FQHE filling factors. It has been 
experimentally established that a,t u = ^ the system is 
metallic. (Here v — (j)Qne/B, where (f>o = hc/e, is the 
mean electron density, and B is externally applied mag- 
netic field.) The seemingly metallic behavior observed 
in transport measurements Q| was confirmed in subse- 
quent surface acoustic wave experiments where it was 
found that contrary to the case of odd-denominator fill- 
ing factors where the conductivity is reduced, at the half- 
filled Landau level the conductivity is in fact, enhanced 
[^. Earlier theoretical attempts to understand the na- 
ture of the V = ^ state, largely pioneered by Halperin ||] 
and later by others [Q-^ remained mostly inconclusive. 
While it was known from those theoretical works that 
the i-statc is compressible, the exact nature of the state 
remained unclear. For example, working with upto ten 
electrons in a periodic rectangular geometry and the ex- 
act diagonalization of the few-electron Hamiltonian in the 
lowest Landau level, Haldane found that the excita- 
tion spectrum is particle number dependent, the ground 
state energy was never at the zero total momentum (con- 
trary to what one expects in a uniform-density liquid), 
and no clear physical picture could be extracted from 
those numerical results. The ground state energy (in fact, 
the lowest energy) was also dependent on the electron 
number and extrapolation of the energies to an infinite 
system led to Eq « -0.465eVe4 |] (here Iq = {hc/eBj 
is the magnetic length). The Laughlin wave function 



(1) 



i<3 



where z = x + iy is the electron position and — 1/m, 
describes a system of particles obeying Bose statistics 
for m = 2 and can not be used for the fermion system 
without any further modification. 

In order to explain the anomalous results at = i, 
a very intriguing theory was proposed by Halperin, Lee, 



and Read (HLR) This theory describes the com- 

pressible even-denominator states in terms of a transfor- 
mation that represents each electron as a Chern-Simons 
fermion carrying even number of fictitious magnetic fiux 
quanta pointing in the direction opposite to the external 
magnetic field. In a mean field approximation (no in- 
terparticle interaction) the average fictitious field cancels 
the real magnetic field and as a result the transformed 
fermions experience no net field. They then form a gap- 
less Fermi liquid. Subsequent experimental observation 
|]lT| of the geometric resonance of the quasiparticle cy- 
clotron orbits with the acoustic waves, and similar ge- 
ometrical resonances found in antidot arrays, indicated 
the existence of a Fermi surface a.t v — ^. These experi- 
ments provided strong support for the theoretical picture 
of HLR. However, fluctuations beyond the mean field the- 
ory, which is expected just to renormalize the Fermi liq- 
uid parameters, instead found to cause divergences and 
the situation has not improved much from there yet [p^ . 

Parallel to the above approach, there is an ongoing 
effort to develop a microscopic approach to v = ^ based 
on the idea of having an improved Laughlin-like wave 
function as a starting point. One way to do that is to 
include the antisymmetrization operator to the Laughlin 
state and have a trial wave function |13,14| 



* = Vlll det M W {z, - Zjf exp ■ 



k 



(2) 



Here Vlll is the lowest Landau level projection opera- 
tor. The matrix M has elements that are plane waves. 



Ikl < kp. Fori/ 



the Fermi wave vector 



is kp = [47rne/s]^ = l/(-\/s£o) where s is the spin degen- 
eracy. For a fully spin-polarized system s = 1. Because 

d 

of the projection operator, Zi — > 2— — , and therefore the 

azi 

plane wave factors act as operators on the Jastrow factor 
where, as a result, the zeroes of -0^ are displaced p3| , p^ . 
The wave function (2) is supposed to have the right statis- 
tics and right correlations to describe the Fermi liquid 
properties at v = i, and is found to provide a good de- 
scription of a small size system at v = ^ p^ . However, in 
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those numeriucal studies of the few-electron systems the 
"Fermi" pair-correlation function on a sphere was found 
to have distinct long-range type oscillations unlike the 
dominant short-range order present in a fluid and also 
not present in the Laughlin ("Boson") state. Further, it 
is reasonable to question the reliability of a few-electron 
system result when we are to describe a gapless Fermi 
liquid. The pair-correlation function and the structure 
function for the state (2) are the most essential build- 
ing blocks for any further development in the theory of 
a compressible fluid. The nature of the correlation func- 
tions in the thermodynamic limit, the effective mass, and 
collective excitations, which are related to those correla- 
tion functions, therefore need careful attention 

In this work, we have attempted to fill in for some 
of those open questions by appealing to the ingenuity 
of the original Laughlin approach where one is able to 
map the electrons onto a classical plasma and make use 
of the established formalism to calculate various physi- 
cal quantities. To develop such a many-body scheme to 
deal with the wave function (2) which includes the pro- 
jection operator is however a non-trivial problem. In our 
approach, we start with a simplified problem and drop 
the projection operator from (2): the wave function is 
then a simple product of the Slater determinant and the 
Laughlin function for the ^-state 



= det M-0L = 



(3) 



Our justification for that somewhat radical step is that, 
in our choice of the trial wave function the only job 
of the Slater determinant $ is to make antisymmet- 
ric. When ^ — 1, the wave function describes the boson 
(Laughlin) fluid and when ip^ = 1 the wave function 
describes the non-interacting Fermi system. Although 
we have no longer any explicit projection to the lowest 
Landau level, the Laughlin wave function, in particular 
the analytic part of the wave function already describes 
the correlations in the lowest Landau level. Also, since 
we are interested primarily in the correlation functions, 
structure functions etc., it should perhaps still be an ac- 
ceptable step to drop the projection operator, especially 
since the form of |<I>p is chosen to be of the same form as 
l^ip (see below). We wish to add here that for a Fermi 
liquid in the absence of a magnetic fleld, a division of 
labor as for the two functions in (3) is entirely justifled. 

Once the choice of the wave function (3) is made the 
next question is, how do we deal with We have al- 
ready stated that we are mostly interested in the pair- 
correlation functions where information about is all 
that need to be known, or more specifically, we need 
to constrain |$p to be positive definite. One available 
choice in the literature |jl^ which is quite successful in 
describing the correlated electron systems in the absence 
of an external magnetic field is to write 



i<3 



(r) = e 



(4) 



where the set of spin coordinates is denoted by a. This 
means that we expand and retain only the two-body 
term which is then approximated by a Jastrow-type func- 
tion. This allows us to write the square of the total wave 
function as 



1^ (ri, • • • , rjv)|' - n e- )+"^ 
and the corresponding pair-correlation function 



(5) 



(6) 



where N is the particle number. 




FIG. 1: The pair-correlation function g[x) as a function 
oi X — r/R^R = ^/2m£Q for the non-interacting system 
[5/ (a;)] [Eq. (7)], the Laughhn state (1) for a Bose system 
{9l{x)], and the "Fermi" state (3) grix). A blowup of 
the region around unity is given as inset. 



The advantage of our choice of (5) is that one can now 
use established methods like the celebrated mapping of 
Laughlin's wave function (1) to a one-component classical 
plasma |^,^ which determines the u(r) or equivalently 
the pair-correlation function g{r). In order to perform 
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similar calculations for Ut{r) = u{r) + ui{r), we first have 
to determine u/(r) and then follow the plasma analogy 
to solve for ut(r). In a completely degenerate ideal (non- 
interacting) two-dimensional Fermi system the exact two- 
body radial distribution function gi{r) can be calculated 
to be ITtI 



gi{r) = l-[2Ji{kFr)/kFrY 



(7) 



where Ji{kpr) is the Bessel function of the first kind 
of order one. Since we are considering a fully spin po- 
larized system, gi{r) vanishes at the origin due to the 
Pauli exclusion principle (Fig. 1). The corresponding 
ideal gas static structure function obtained from the two- 
dimensional Fourier transform of (7) is jl^ 



Si{n) 



1, 



K (l 



K <1 
K > 1 



(8) 



where k — k/2kp. For small k, Si{k) increases linearly 
with K (Fig. 2). 
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FIG. 2: The static structure factor S{q) as a function of 
q = kR, for the non- interacting system [^/(g)] [Eq. (8)], 
the Laughhn state (1) for a Bose system {SL{q)], and the 
"Fermi" state (3) Spiq)- A blowup of the region around 
unity is given as inset. 



Given these two functions, u/(r) now can be obtained 
from inverting the hypcrnetted-chain (HNC) equations 



ui{r) = - In .9,(0 + <^ ^ Vrt.^ 

He Sl(k) 



(9) 



F{k) = 27r / rdrF{r)Jo{kr) 
Jo 

F{r) = ^ f kdkF{k)Jo{kr), 



where jQ{kr) is the Bessel function of the first kind of 
order zero. 

The implication of Eqs. (7)-(9) is that in the two-body 
level, (j>{r) in Eq. (4) can be chosen such that we repro- 
duce the exact two-body radial distribution function and 
the static structure factor corresponding to the full de- 
terminant. We note that ui{r) is a numerically strictly 
decreasing function, but it is a long-ranged function (be- 
cause for small q, ui{q) = J-[uj{r)] ^ q^^)- There- 
fore one needs to take special care about the long- and 
short-range behavior of ut (r) ||l^] while solving the HNC 
equations for the one-component plasma. In our numeri- 
cal calculations we have used the dimensionless variables 
X = r/R and q = kR where R = V2m£o is the ion-disk 
radius |^ . The method of deriving the HNC equations for 
Ut{x) is similar to that for u{x) and with proper choice of 
the long- and short-range functions [§,0 a numerically 
rapidly convergent set of equations are obtained which 
lead to g{x) and its Fourier transform, S{q). In the case 
of ui{x) ~ 0, the pair-correlation function gL^x) for the 
Laughlin state is plotted in Fig. 1. The ground state 
energy coresponding to that state is El = — 0.480e^/e£o- 
The pair-correlation function corresponding to state (3) 
is denoted by gwix) in Fig. 1. Clearly, the "Fermi hole" 
is not much affected by the introduction of the Laugh- 
lin correlation function, but marked deviations of gpix) 
from gi{x) occur near the maximum of gL^x). 

The static structure functions S{q) vs q for the various 
cases are shown in Fig. 2. As with the pair-correlation 
functions, Si{q) corresponds to the ideal system result, 
Spiq) is the present result for the Fermi-fluid state (3) 
at J/ = i and S'l('?) is the structure function for the 
Laughlin state (1). For small q, it has been predicted 
that S{q) cx in the Fermi fluid a,t ly = ^ Q. Simi- 
lar behavior is observed here for Sp{q). The difference 
between the present results and the Laughlin results, 
gdiff = gpix) - gL{x) and Sdiff = Spiq) - S'l(9), are 
plotted in Fig. 3. These are oscillatory functions with a 
rapidly decreasing amplitude. In finite-size system calcu- 
lations [0, a sinusoidal oscillation in gdiffix) was taken 
as an indication of a Fermi fluid behavior. We note that 
the correlation functions for the non-ideal system show 
much less oscillations around unity in accordance with 
the properties of a uniform density fluid (inset of Fig. 
1) and therefore the difference in correlation functions 
is also rapidly damped. Interestingly, Sdiff (k) develops 
a positive peak slightly below k ^ 2kp and a negative 
peak beyond that k. Finally, we find the ground state 
energy for the state (3) and for the Coulomb potential 



where J- denotes the two-dimensional Fourier transform to be Ep — —0.448 e^/e^Oj which is very different from 
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the energy of the LaughUn state El but very close to the 
energy value Eq, extrapolated for an infinite system from 
the finite-size system results mentioned in the introduc- 
tion. This agreement between the energy of the state (3) 
and the estimate Eq is a strong indication that our Fermi- 
liquid description has the right correlations and correct 
statistics needed to describe a Fermi liquid behavior at 
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FIG. 3: The difference between the "Fermi" and "Bose" 
functions (a) gdiff(x) vs x and (b) Sdiff{q) vs q. 



In summary, our simplified choice (3) of the ground 
state wave function for the Fermi fiuid state at v = ^ 
has led to a microscopic approach where we can calcu- 
late the physical quantities in the thermodynamic limit. 
The pair-correlation function, the structure function, and 
the ground state energy are in good agreement with the 
expected behavior at v — ^. This approach can also 
be suitably modified to calculate the one-body density 
matrix and the nature of the off-diagonal long-range or- 
der (ODLRO) which will provide more information 
about the Fermi nature of the proposed state. In de- 
fense of our choice of (3), we should mention that the 
wave function which would result from the operation of 
"Plll will be a wave function in the lowest Landau level, 
and therefore be similar to the Laughlin-like wave func- 
tion (but with correct statistics). Hence our choice of 
(5), which is formally similar to the Laughlin approach 
should be a suitable approximation for the full wave func- 
tion (2). This is supported by our numerical results pre- 
sented here. 
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